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Abstract 

We discuss a novel mechanism for J /if/ production in nuclear collisions arising due to the high 
density of gluons. We demonstrate that gluon saturation in the colliding nuclei is a dominant 
source of J /if/ suppression and can explain its experimentally observed rapidity and centrality 
dependence. 



1 1. Introduction 

2 The mechanism of J/ifi production in high energy nuclear collisions is different from that in 

3 hadron-hadron collisions flll2l. Consider first the J/ifi production in hadron-hadron collisions. 

4 The leading contribution is given by the two-gluon fusion (i) G + G — > J I if/ + soft gluon, see 

5 Fig. [i]A. This process is of the order 0(a 5 s ). The three-gluon fusion (ii) G + G + G — > J/if/, 
e see Fig. [l}B, is parametrically suppressed as it is proportional to 0(a 6 s ). However, in hadron- 
7 nucleus collisions an additional gluon can be attached to the nucleus. This brings in an additional 
s enhancement by a factor ~ A 1 ^ 3 . If the collision energy is high enough, the coherence length 

9 becomes much larger than the size of the interaction region. In this case all A nucleons interact 

10 coherently as a quasi-classical field [3] HI El |6) . In the quasi-classical approximation a 2 A 1/3 ~ 1. 

11 Therefore, the three-gluon fusion is actually enhanced by l/a s as compared to the two-gluon 

12 fusion process. Similar conclusion holds for heavy-ion collisions (we do not consider any final 

13 state processes leading to a possible formation of the Quark-Gluon Plasma). This approach has 

14 been previously applied in Ref. (7) to 7/^ production in dAu collisions, although the nuclear 

15 geometry was oversimplified. We note a reasonable agreement of this earlier approach with the 

16 dA data. 

7 2. New mechanism of J /iff production 

Fig.[T]-B represents the contribution of the order (a 2 A 1 ^ 3 ) 2 . In general, there must be an odd 
number of gluons connected to the charm fermion line because the quantum numbers of Jfi/i and 
of gluon are 1 . Therefore, each inelastic interaction of the cc pair must involve two nucleons 
and hence is of the order (a 2 A 1 ^ 3 ) 2 ", where n = 1,2, . .. , A/2 is the number of nucleon pairs. 
To take this into account it is convenient to write the cross section as the sum over all inelastic 
processes (labeled by the index «). This sum involves only even number of interactions. For a 
heavy nucleus A » 1 and for N c » 1 the calculation can be significantly simplified. The final 
result can be written in terms of the saturation scale Q 2 A ||8] [10) . Its value was extracted from 
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Figure 1: The process of inclusive production in hadron-hadron (Fig.[TjA) and in hadron-nucleus collisions (Fig.QJ 
B). 

the fit of the multiplicities of nuclear reactions at RHIC ATI [121 and from fits of the Fi structure 
function in DIS 02] [H E] [BD . The main result reads fl]|2l 



C d _E^l J d 2 s ^ (£) Ta _ ^ _ £ ^ ^ > £ + Q 2^ ^ b _ _ ^ 

X°° r fi 

d^t? K 2 (0 exp[-^ {QIaMusJ + Q 2 sA2 (x 2 ,b- s)) , 



(1) 



where JC1/2 = ^= e* 7 . Eq. jli is derived in the quasi-classical approximation which takes into 
account multiple scattering of the cc pair in the cold nuclear medium. At forward rapidities at 
RHIC and at LHC the gluon distribution function evolves according to the evolution equations 
of the color glass condensate. Inclusion of this evolution in the case of J /iff production presents 
a formidable technical challenge. Therefore we adopt a phenomenological approach of lfl2l [T71 
in which the quantum evolution is encoded in the energy/rapidity dependence of the saturation 
scale. 

In our numerical calculations we take explicit account of the impact parameter dependence 
of the saturation scales of each nucleus. We employ the Glauber approximation and assume that 
the nucleon radius is much smaller than the nucleus radius. The overall normalization constant 
C includes the color and the geometric factors C 2 F /{An 2 a s S p ) where S p is interaction area in 
proton-proton collisions. C also includes the amplitude of quark-antiquark annihilation into 
J/t// and a soft gluon in the case of pp collisions. This amplitude as well as the mechanism of 
Fig.[T]-A have a significant theoretical uncertainty. Therefore, we decided to parameterize these 
contributions by an overall normalization constant in ([TJ. 

The rapidity distribution of J/t/z's in pp collisions, the factor dN pp /dY appearing in Eq. [T] 
is fitted to the experimental data given in [18] with a single gaussian. In figure |2|a) the results 
provided by Eq. ([T]l are then compared to the experimental data from PHENIX Collaboration 
|[T9l for Au-Au collisions at = 200 GeV. The global normalization factor C is found from 
the overall fit. There are no other free parameters. The agreement of the theoretical results 
with experimental data is reasonable. It is evident that the effect of the gluon saturation on the 
J/t// rapidity distribution in nucleus-nucleus collisions is to make its width a decreasing function 
of centrality. The distribution in the most central bin is significantly more narrow than in the 
peripheral bin. 

It is important that we describe well the data in the semi-peripheral region. This ensures that 
our model gives a good description of the J /iff production in dAu collisions. We also note that an 

2 




earlier approach [7| in which the same model was employed (albeit with an oversimplified nu- 
clear geometry) provided a reasonable description of the data. Still a more detailed investigation 
is required which takes into account the exact deuteron and gold nuclear distributions. This will 
allow a model-independent fixing of the overall normalization constant C. We plan to present 
such an analysis in the near future. 

To emphasize the nuclear dependence of the inclusive cross sections it is convenient to in- 
troduce the nuclear modification factor defined in a usual way. In Fig.[2|b) we plot the result of 
our calculation. The nuclear modification factor exhibits the following two important features: 
(i) unlike the open charm production, J /if/ is suppressed even at y — 0; (ii) cold nuclear matter 
effects account for a significant part of the "anomalous" J /if/ suppression in heavy-ion collisions 
both at y = and y = 2. 



3. Conclusions 

The main results of this paper are exhibited in Fig. [2] It is seen that the rapidity and centrality 
dependence of J /if/ production are reproduced with a reasonable accuracy even without taking 
into account any hot nuclear medium effects. This observation allows to conclude that a fair 
amount (and perhaps most) of the J/if/ suppression in high energy heavy-ion collisions arises 
from the cold nuclear matter effects. In other words, J /if/ is expected to be strongly suppressed 
even if there were no hot nuclear matter produced. 

The reason for J/i/j suppression at mid-rapidities is that the multiple scattering of cc in 
the cold nuclear medium increases the relative momentum between the quark and antiquark, 
which makes the bound state formation less probable. It was proven in |20| that unless quantum 
log(l/x) corrections become important, the inclusive gluon production satisfies the sum rule that 
requires the nuclear modification factor to be of order unity. Similar sum rule holds for heavy 
quark production but fails in the case of a bound states, such as J/i/t. 

We realize that although our calculation gives the parametrically leading result at high gluon 
density, other production channels involving the gluon radiation in the final state and the color 
octet mechanism of J /if/ production may give phenomenologically significant contributions. 
These are likely to become the leading mechanisms in the peripheral collisions where the strength 



of the gluon fields is significantly diminished. However, we believe that our main results are ro- 
bust for central high-energy collisions of heavy ions. 
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